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Abstract 

We present a theory of the finite temperature thermo-electric response functions of graphene in 
the hydrodynamic regime where electron-electron collisions dominate the scattering. In moderate 
magnetic fields, the Dirac particles undergo a collective cyclotron motion with a temperature- 
dependent relativistic cyclotron frequency proportional to the net charge density of the Dirac 
plasma. In contrast to the undamped cyclotron pole in Galilean-invariant systems (Kohn's theo- 
rem) , here there is a finite damping induced by collisions between the counter-propagating particles 
and holes. This cyclotron motion shows up as a damped pole in the frequency dependent conduc- 
tivities, and should be readily detectable in microwave measurements at room temperature. We 
also compute the large Nernst signal in the hydrodynamic regime which is significantly bigger than 
in ordinary metals. 
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I. INTRODUCTION 



In the absence of gate voltage or external impurities, graphene is a quantum critical sys- 
tem , whose quasiparticles are governed by a relativistic massless Dirac equation 
in two spatial dimensions 0, Q, Isl- (By 'quantum critical' we mean here a system in which 
the inelastic scattering rate from electron-electron interactions is of order kBT/h (up to 
logarithms), and not models of localization transitions of free electrons Experimental 
realizations of relativistically-invariant systems are rare, and so it is of interest to study phys- 
ical phenomena which rely on the Dirac nature of the electrons. In this paper we focus our 
attention on the possibility of observing a relativistic cyclotron resonance in the (collective) 
electronic motion in graphene. This resonance shows up as a bump in the frequency depen- 
dence of all thermoelectric response functions. As we will discuss in detail, the cyclotron 
resonance frequency has clear signatures in its dependence upon field (B), temperature (T), 
and electron density (p) which distinguish it from the cyclotron resonance of non-relativistic 
electrons. Moreover, the relativistic cyclotron resonance is intrinsically damped by electron- 
electron interactions: this damping arises from collisions between electrons and holes which 
execute cyclotron orbits in opposite directions. In contrast, the damping of the cyclotron res- 
onance in metals arises primarily from impurity scattering: Kohn's theorem 10| implies that 
electron-electron interactions do not broaden the cyclotron resonance in Galilean-invariant 
systems. This theorem applies to metals with a single parabolic (non-relativistic) band, 
which is a reasonable approximate description for many simple Fermi surfaces. However, it 
fails for semimetals such as graphene both due to the linear, relativistic dispersion and the 
presence of two bands. 

We shall be interested here in the cyclotron resonance in a hydrodynamic, collision- 
dominated regime, where disorder plays only a minor role. This regime exists at high 
temperature and is defined by the requirement that the vast majority of collisions arise 
from electron-electron interactions. This assures local equilibration before scattering from 
impurities occurs. The second requirement for hydrodynamics to apply is the following: 
The rate of deflection of linearly propagating, thermal quasiparticles due to the magnetic 
field must be small compared to the inelastic scattering rate, t-~1^, arising from interactions. 
The same restriction must hold for the frequency of the external driving fields, cu <^ t^"^^. 
The quantum-critical nature of graphene [l| implies that the mean time between collisions 
is of order /i/fc^T. Since we require that this inelastic collision time be shorter than the 
elastic scattering time from impurities, we cannot allow T to become too small, and will find 
that room temperature is suitable for observing the physics we are interested in. Further 
conditions controlling the range of parameters are discussed below in Section IIIII 

The cyclotron resonance is formally defined as the pole of the thermo-electric response 
functions closest to the real axis in the complex frequency plane. The expressions for the 
cyclotron frequency (its real part), Uc, and the damping frequency (its imaginary part), 7, 
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are the same as those in Refs. [ill, Il2l : 



Here vf is the 'velocity of hght' for the Dirac fermions, which was experimentally measured 
to be jisl, Q, ~ 1-1 X 10^ cm/s ~ c/300. The density of electrons p is defined so that 

p = for undoped graphene. e and P are the thermodynamic energy density and pressure of 
the Dirac plasma, respectively, which are also measured with respect to undoped graphene. 
In the relativistic regime where the temperature exceeds the chemical potential the 
energy density and the pressure grow with the third power of temperature and assume 
typical values of e, P ~ {kBTf/ivphf = 3.28 x lO^^ ^ TlggmeV/cm^ where T300 = T/300K. 
The coefficient ag will be seen to cinSG clS cL transport parameter in the hydrodynamic 
description of a relativistic fluid which cannot be fully determined by thermodynamics and 
hydrodynamics alone. It can, however, be computed in a microscopic approach, as recently 



carried out in Ref, 15(]. The parameter ctq has the units of an electrical conductivity and, 
as will become clear from later formulae, it describes the part of the d.c. conductivity which 
is independent of impurities, deriving solely from interactions, c.f., Eq. (15.11) . At particle- 
hole symmetry, ctq coincides with the finite conductivity of a clean system in zero field, 
B = 0. In the relativistic, collision-dominated regime which is of foremost interest here. 



one has ctq ~ [16|] where a = jj— ~ is the fine structure constant of graphene 



characterizing the strength of Coulomb interactions, and is the dielectric constant due 
to the adjacent media. In general, ctq it is a scaling function of p/T. The dependence of 
the thermodynamic variables and the transport coefficient aq on temperature and chemical 
potential are further discussed in Section IIII A[ 

Physically, the cyclotron resonance is due to the tendency of fiuid elements of the electron- 
hole plasma to undergo a circular motion at frequency Uc- This frequency results as an 
average over the left- and right-circulating orbital motions of thermally excited electrons 
and holes, which collide with each other many times before they would be able to complete 
a (non- interacting) cyclotron orbit. For this reason, Uc is proportional to the excess charge 
density p, and vanishes at the particle- hole symmetric Dirac point where p = 0. Nevertheless, 
a bump in the frequency dependent response around 00 = survives also in this case. Note 
that it broadens rather rapidly with decreasing T as 7 ~ T~^. 

In the non-relativistic regime \p\ ^ T, where graphene turns into an ordinary Fermi liq- 
uid, the cyclotron resonance tends to the semiclassical value Uc = eVpB/{cp) corresponding 
to the cyclotron frequency of non-interacting Dirac fermions at the chemical potential p. 
This should indeed be expected since all thermally excited quasiparticles share essentially 
the same non-interacting cyclotron frequency which is not altered by the interactions. 

From microscopic transport theory one finds that aq decreases as {T / pY in the non- 
relativistic Fermi liquid regime [l5] . Accordingly, the intrinsic damping 7 decreases and the 
cyclotron resonance becomes increasingly sharp for small T / p. We note that in addition 
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to the intrinsic damping represented by 7, there will also be extrinsic damping from im- 
purity scattering, as is discussed in Section lllli The complete expression for the frequency 
dependence of the conductivity across the cyclotron resonance is given in Eq. fl5.6p . and is 
illustrated in Fig. [2l 



Recent experiments 13|, ll8| have observed a "non-hydrodynamic" cyclotron resonance, 
in a regime of strong magnetic fields in which the Landau levels can be resolved. As we will 
see below, we are discussing here the different regime of weak fields and high temperatures, 
to which we hope the experiments will be extended. 



II. MODEL OF GRAPHENE 

We consider a single sheet of graphene described by the Hamiltonian 

H = Ho + H, + H^,,, (2.1) 
where the low energy tight binding part is 

Ho = -J2 Yl I d^\^\y{iVFB-y + ii^'^sy] , (2.2) 

s V=K,K' 

with the Fermi velocity vp^ the two component wavefunction sy describing the amplitude 
on the two sublattices for electrons with spin s and momenta close to one of the two Fermi 
points V = K' . A magnetic field is introduced as usual via minimal coupling. 
The 1/r Coulomb interactions take the form 



^1 = ^EE / d^d^^ly{-K)^>,y{^) 
s'V 



The term ifjis describes the presence of weak disorder which induces elastic scattering and 
thus weak momentum relaxation at a rate r proportional to the impurity concentration. In 
the high temperature regime which we are focusing on, the elastic scattering rate is smaller 
than the inelastic scattering rate and can be taken into account as a perturbation. 



III. THE OBSERVABILITY OF THE COLLECTIVE CYCLOTRON RESONANCE 
IN GRAPHENE 

As discussed in Ref. the applicability of hydrodynamics to the system (12. ip requires 
that the magnetic field B be weak, so that the Landau-quantization of thermal excitations 
is not discernible, i.e.. 



E, = hvF\l ^ keT, (3.1) 
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where Ei is the first Landau level for graphene at the Dirac point 

0, 0, HQ- This 

amounts to requiring that 

B « B\T) ^ YeJ^ = X T^o) Tesla. (3.2) 
For the hydrodynamic analysis to hold, one needs a slightly more stringent condition: 

El < ^^^ei ~ «'^bT, (3.3) 

where the last estimate applies to the relativistic regime T < This condition expresses 
the requirement that the inelastic scattering rate, and thus the equilibration rate, should 
dominate the rate by which the magnetic field defiects electrons from their linear motion. 
In the following we will assume a to be of order unity, but we nevertheless will retain a in 
most scaling estimates. 

It will also be useful to express our electron densities in terms of the characteristic density 
of thermal excitations, 



Ah-(|^)'=(l.27xlO-xT3^,o) 



cm-^ (3.4) 



This density should be compared to disorder induced density variations which are of order 



^Pdis ~ lO^^cm"^ ij] varying on typical length scales ~ 30nm. To remain close to quantum 
criticality, ensuring a universal conductivity due to thermal pair creation and annihilation 
processes, the regime (5pdis < Pth is preferred. This suggests measurements at room temper- 
ature or above. 

As mentioned above, the cyclotron resonance occurs at 

eBpvl _ 1 keT B p 
c{e + P) ~ 2$,+p h B*{T)'^ 

1.96xl0^^s~\ (3.5) 



$,+P B*{T) Pth 



where we have used a free-electron equation of state to determine the value of e + P, with 
$£_i_P = {e + P){fivYT~^ being a dimensionless number 0(1), as given in Eq. fl3.13p below; 
we will also comment in Section IIII Al on the effect of interactions on the equation of state. 

The collective cyclotron frequency lies within the hydrodynamic frequency regime if 
hujc ^ ^''^nci ~ cP'kBT where the latter estimate applies to the relativistic, quantum-critical 
regime. This requires 

hjjn 1 B p , , 

<1. (3.6) 



a^ksT 2a^^,+p B*{T) p,^ 

For room temperature and values of B and p as suggested below, this falls into the range of 
microwave frequencies. Thus the resonance should be readily detectable by measuring the 
real and imaginary part of (Txx{^) of a graphene sample in a cavity. 
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For inelastic scattering to be the dominant relaxation process, i.e., one certainly needs the 
impurity scattering time to satisfy r ^ Tinei. However, to observe the relativistic cyclotron 
resonance, one has to require the more stringent condition 

r>lM= - — = — 5x10 s, (3.7) 

\nWc/ isOO ^300 B p 

to ensure that the disorder- induced broadening does not wash out the resonance. Such long 
scattering times can indeed be achieved in high mobility graphene where r ~ 10~^^s is a 
typical value. 0, EH 

Apart from a small impurity scattering one also needs the intrinsic broadening of the 
cyclotron resonance due to electron-electron interactions to be smaller than Uc- 

7 _ B Pth . . 

where $^ = ag/ieyh) = 0(1). 

Note that ( 13. 6p and (13. Sp can be satisfied simultaneously, e.g., with a choice of magnetic 
field B/B*{T) ~ 0(0. 1) and charge density p/pth ~ 0{l). 



A. Scaling functions for thermodynamic variables 

For the evaluation of thermodynamic state variables magnetic field effects can be safely 
neglected since both Zeeman energy and Landau level splitting are significantly smaller 
than T. Further we will also neglect interactions and use the free theory to evaluate ther- 
modynamic quantities. This can be justified by the fact that the Coulomb interactions 
are marginally irrelevant 0, 0] . However, as emphasized by Son [1] , the bare value of the 
interactions are not small (a = 0(1), and there may well be a significant regime of inter- 
mediate scales where the interactions remain significant, and the theory is characterized by 
a dynamic critical exponent 2; 7^ 1. Over this intermediate regime, we have e + P ~ T^+^. 
Below, we neglect such effects, and simply use the free theory to obtain a numerical estimate 
of parameters. 

One obtains straightforwardly the scaling functions 

£ + P = 7^^e+p{i^ = /i/T) = Tpth$.+p(p = p/T), (3.9) 



{hvy 



with 



$e+p(p) =4^^ 



3 f^kdk 



27r 



k k 



ok-jl _j_ gfc+A _|_ \ 



9C(3) 6 log(2)u2 u4 
sv ; ^ h\ ;p + ^ + 0{Pi^). (3.10) 



TT TT 8 TT 
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The prefactor 3/2 derives from the relation P = e/2, and the factor 4 accounts for spin and 
valley degeneracy. The electron density is 



y2 



P 



{hv) 



-,^p{p = P/T) = pth^pip = p/T), 



(3.11) 



with 



kdk [ 1 
4 log(2) fl fi^ fl^ 



1 



+ 



6 TT 240 TT 



(3.12) 



and thus 



$e+p(p = p/pth) 



9C(3) , 37rp 



;;2 



37r3p4 



71 81og(2) 2048 log^ (2) 
3.444 + 1.700p2 - 0.1968p^ + 0{p^). 



(3.13) 



A plot of this function is shown in Fig. [T] We note that it has been computed in the weak- 

25 
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P/pth 

FIG. 1: The dimensionless function $e+p of the density p/pth for non-interacting Dirac fermions. 

coupling limit; if a = 0{1) there will be order unity corrections from the Coulomb interac- 
tions which can be determined from the theory in Ref. Q. 

For the results given in the next section the scaling function for the entropy density 
s = {e + P — pp)/T will also be useful: 

— = $,(A) = <l>,+p(p) - p<l>,(p). (3.14) 
Pth 
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2 

Finally we mention that the transport coefficient ctq = ^$o-(At/T) is itself a scaling 
function of the ratio between chemical potential and temperature. In the relativistic regime, 
T > /i, it can roughly approximated by a constant, whereas in the non-relativistic, Fermi 
liquid regime, T < /i, a microscopic study shows that ctq decreases as (T/yu)^. The precise 
scaling function has been obtained in Ref. 
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IV. DERIVATION: HYDRODYNAMICS OF A RELATIVISTIC FLUID WITH 
COULOMB INTERACTIONS 



Here we discuss the magnetohydrodynamics of a relativistic fluid HI] in the presence 
of Coulomb interactions. Our general results for response functions in the hydrodynamic 
regime are then readily applied to graphene close to its Dirac point. The characteristic 
velocity vp ^ c determines the relativistic dispersion, and the charge of (anti-)particles in 
the fluid is ±e. In the following we employ units with vp = e = h = \. 

Because the Coulomb interactions spread with the speed of light c'^ vp, they can be 
considered instantaneous. They obviously break the relativistic invariance of the fluid by 
singling out the lab frame as a particular reference frame, which we will eventually work 
below. 

The stress energy tensor, T^^ ^ and the current, J'^, of a relativistic fluid are expressed 
in terms of the three- velocity of the fluid element, = {l,v)/\^l — v'^ (u^u,,, = —1, with 
= (1, 0, 0) in the local rest frame for which the energy flow vanishes): [lll.l22| 

T^-" = (^e + P)u^'u'' + Pg^" + T^'\ (4.1) 
Ji" = pu^ + u^". (4.2) 

These expressions include dissipative terms u'^, t^^ which account for heat currents and 
viscous forces, respectively, and will be derived below. Without the viscous terms, the stress 
energy tensor is a diagonal matrix in the rest frame of the fluid element with the pressure in 
the space-like entries, and energy density in the time-like entry. In the lab frame (to lowest 
order in f*) the components read 

T°° = £, (4.3) 
= = (e + P)v\ (4.4) 

rpij ^ P5'i + T'\ (4.5) 
J° = p, (4.6) 
r = pv' + v\ (4.7) 

Here e, P, p, are functions of the local chemical potential, /i(r), the local temperature T(r) 
and the magnetic field -B, as given in the previous section for graphene. Thereby the static 
Coulomb potential created by an inhomogeneous charge distribution is incorporated in /i(r). 
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The conservation laws for charge, energy and momentum read: 



d,J^ = 0, 



7' 



(4.8) 
(4.9) 



where the electromagnetic field tensor, 



( Ey\ 

-E^ B 
V -Ey -B j 



(4.10) 



contains a self-generated, spatially varying electric field due to the inhomogeneous charge 
density of the system itself: 



E = -V0, 

d^yU{x-y) [p{y) - p] 



[X 



(4.11) 



with 



f/(r) = - 



r J (27r)2 
E{k) = -ikUkp{k), 



Uk exp[iA; ■ r] 



2'Ka 



(4.12) 



where the uniform background charge density p has been subtracted in (14. lip . 
In the lab frame the linearized conservation laws read more explicitly: 

+ V ■ J = 0, 

dtS + V ■ [{e + P)v] = pv-E, 
{dt + r-') {e + P)v + VP = BeJ + pE, 



(4.13) 
(4.14) 
(4.15) 



where in addition we have included a momentum relaxation time r due to weak impurity 
scattering. Note that the latter also breaks the relativistic invariance, and accordingly the 
relativistic hydrodynamics should only be expected to hold as long as r^^ <^ t^I^. The 
above set of equations is closed by the constitutive equation for the current 



J = pv + u. 



(4.16) 



A. Heat current and viscous terms in a magnetic field 

With the help of the thermodynamic relations 

P + e = sT + fip, de = Tds + fidp, 



(4.17) 
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and employing the conditions [22| m^i/'^ = u^t 



0, the energy conservation law 



{Ufj_dyT^^ = u^F'^'^Jij) can be rewritten in the form of a law of entropy production, 



It is natural to interpret the left hand side as the divergence of an entropy current, 



T 



(4.18) 



(4.19) 



which, by the second law of thermodynamics, has to be positive. For small velocity deriva- 
tives d^u^ and electromagnetic fields B, E, the dissipative terms should be linear in these 
perturbations, and hence one deduces the form 



(4.20) 
(4.21) 



where ag is a conductivity (of order e'^/h), and rj and ( are the shear and bulk viscosities, as 
will be clear from the expression (14.231) given below. For graphene at the Dirac point, V, C y 
T'^ jo?. These can be computed by solving a linearized kinetic Boltzmann equation 16l. |23|. 
The relativistic expression (14.201) reduces to a spatial vector in the lab frame 



(4.22) 



V = -cxq [Vfi - ^jj - Bev-Ej , 
while the viscous contribution to the stress energy tensor takes the familiar form 



r ■ 



(4.23) 



B. Linear thermoelectric response 

The thermoelectric transport coefficients describing the current (J) and heat current (Q) 
response to electric fields and temperature gradients are defined by the relation 




(4.24) 



where a, d, a and k are 2x2 matrices acting on the spatial indices x, y. Rotational invariance 
in the plane imposes the form 

a = a^^i + (T^yi, (4.25) 

where 1 is the identity, and e is the antisymmetric tensor i^y = —^yx = 1- cFxx and a^y 
describe the longitudinal and Hall conductivity, respectively. An analogous form holds for 
the thermoelectric conductivities d, a which determine the Peltier, Seebeck, and Nernst 
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effects, as well as for the matrix k which governs thermal transport in the absence of electric 
fields. The latter applies to samples connected to conducting leads, allowing for a stationary 
current flow. In contrast, the thermal conductivity, /t, is defined as the heat current response 
to — VT in the absence of an electric current (electrically isolated boundaries), and is given 
by 



K 



K — Taa a. 



(4.26) 



C. Response functions from hydrodynamics 



We will now use the conservation laws ( ]4.13|l4.14|l4.15p to solve for the slow relaxation 
dynamics towards equilibrium, starting from a small initial long-wavelength perturbation. 
From the full solution of the relaxation dynamics in linear response approximation one can 
then determine the thermo-electric response functions in the hydrodynamic regime. 

As dynamic variables we choose T, /i (which includes the static Coulomb potential), and 

and v^, and write 



/i(r, t) = /i + 6fi{r,t), 
T{r,t) = T + 6T{r,t). 



(4.27) 



The other variables, e, P, and p are constrained by local thermodynamic equilibrium to 
have the form 



p{r,t) = p + 6p = p + 



dp 
dp 



5p + 



T 



dp 
df 



e{r,t) 
P{r,t) = 



= e + Se 
P + 5P 



de_ 
dp 



T 



ST, 



ST, 



P + pSp + sST, 



(4.28) 



where Sp is the deviation of the total electrochemical potential form the equilibrium value 
p, and dp/dp\T = X is the susceptibility. Note, that even though for graphene at the Dirac 
point {p = 0) x{k) ^ k at T = 0, the susceptibility is always finite, x ^ T, for the 
hydrodynamic regime considered here (T > 0). 



Following the technique of Kadanoff-Martin 2j] to derive hydrodynamic response func- 
tions, we prepare the system in a state of local equilibrium as characterized by slowly varying 
initial conditions T{r) and p{r), as created by an external electric field. Local equilibrium 
implies that the screening of an inhomogeneous charge density must be selfconsistently built 
into the initial conditions. The relation between the actual initial variation of the chemical 
potential, Sp^, and both an applied external field E'^^^ = —V {Sp'^^^) and temperature de- 
viations ST^, is therefore not entirely trivial: The total chemical potential p{r) is the sum 
of the externally applied potential Sp^^^, and the Coulomb potential created by the total 
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induced charge density. In /c-space, the Thomas-Fermi (RPA) screened Coulomb potential 
induced by the two perturbations leads to the chemical potential 



1 + Uk dp/dfi\ 




T 



ST\k) 



(4.29) 



After a Fourier transform in space and a Laplace transform in time, Eqs. f l4.13tl^rT5|) 
together with ( 14.121) take the form 




5fx + 



T 




5fx + 



de_ 
df 

dp 
df 

dp 
dp 



6T -k{e + P)v\\ 



5p' + 



T 



de_ 
df 



6T° 



(4.30) 



ST^ - k {pv\\ + GQ [Bv^ + Ek]) + lOQe (6p - ^6T 



T 



dp 
df 



6T^ 



where 



u + + P)v\\ - k{p5p + s5T) - ipBv± + ioqE'^vw - ipEk + ik'^{'n + C)^l| = + ^)^|| , 

^ j (e + P)v± + kaqB [dp - ^5T^ + ipBv\\ + igqB [Bv^ + Eu] + ik'^r^v^ = i{e + P)4 



Ek = -ikUk 



5p + 



T 



dp 
df 



6T 



(4.31) 



In these expressions, 6p^, 6T^, 5f y and 6v^ are the initial values (depending on the wavevec- 
tor k), while 6p, 6T, 6v\\ and 6v± are functions of k and u. The projections of v parallel and 
orthogonal to k are v\\ = v ■ k/k, y \ = k/k ■ ev. 



As shown in detail in Refs. [ill . |2j] , the linear response functions can be obtained in full 
generality and in closed form from the solution of the above equations, e.g., a^xi^^jk) = 
ujJ\\{uj,k)/5p^, where J\\ = J ■ k/k. 

We have restricted ourselves to the lowest non-trivial order in an expansion in k/uj, 
where Coulomb effects become visible. This will be sufficient to describe the response to 
external perturbations of electromagnetic origin (e.g., microwaves) for which one always has 
ck/uo < Y^e^. When applying the relativistic hydrodynamics to graphene with a characteris- 
tic velocity f/c ~ 1/300, this implies the relation vk/oj < {v/c)^/e^ <^ 1, which justifies the 
expansion to lowest order in k/uo. The full fc, cj-dependence can be obtained in closed form. 
However, it is very involved and does not contain much more physical information, so we do 
not report it explicitly here. 
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V. RESULTS FOR THE THERMO-ELECTRIC RESPONSE 



A. Limit of vanishing magnetic field 



In the limit of vanishing magnetic field -B ^ 0, the transverse response vanishes for 
symmetry reasons. The longitudinal transport coefficients take the relatively simple forms 



a. 



xx\^t k', B 



iuj, k: B 



i^xxis^i k] B 



0) 
0) 

0) 
0) 



a. 



P + el — iuoT _ 
uj,k;B = 0) + O{k'^ 

fi Sp T 



ik 



1 - 27ra— \(JQ + - 
CO \ 1 



iujT P + e 



"^«T + P + el 



lUJT 



1 - 27ra- 



ik 



UJ 



T 



P' 



P' 



T 



T ^ P + el-iujTj I UJ 

UJ i 1 — lUJT 

(P + ef 



ik 

2na— aq + 



iujT P + e 

■,2 



+ 0(fc^),(5.1) 
(5.2) 



p 



1 — iujT P + e 



1 



T p2 + [aQ/T){P + e){l 



IUJT) 



+ 0{k' 



(5.3) 
(5.4) 



Note that all response functions contain a piece proportional to cxq which is independent 
of the impurity scattering time r and thus is solely governed by the universal Coulomb 
interactions. A second term proportional to t/(1 — iur) has the form of a classical Drude- 
like term which is due to the slow relaxation of the "momentum mode" 16||. This is an 
excitation of the electron-hole liquid which cannot relax due to Coulomb interactions because 
of their translational invariance. As one may expect for weak impurity concentration, the 
term contributed by the momentum mode grows linearly with the impurity scattering time 
r. 

Note the simple structure of the leading k dependence of the response functions: It only 
depends on Coulomb interactions (via a), while the viscosities ri,( do not appear at this 
order. Including finite k introduces a simple factor 



27ra- 



ik 



UJ 



P' 



1 — iujT P + e 



27ra—a^^{uj, k = 0), 



(5.5) 



which multiplies the A; — result. Interestingly, the thermal conductivity k, (in the absence 
of currents) is not affected by Coulomb interactions to lowest order. 
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B. Magneto-transport 



Most interesting for our purpose is the response in a weak magnetic field for which we 
obtain: 



{uo + i/r) {uo + i/r + 27 + iooHi) 

' 2 

(co' + i/r + i7) - uSl 



(5.6) 



1 — 27ra — oq 



ik (uj + i/r) {uj + i/r + i'j + iuj'^/'y) 



LU 



{uj + i/T + 27)^ - ^ 



+ 0{k'), 



a^yiuj, k) 



p ujI + 7^ + 27(l/r - iuj) 

X 



(5.7) 



ik {uj + %l r) (cj + z/r + Z7 + luJlj 7) 
i — ZTTO; — an — 
a; 



+ 0(fc'). 



(cj + z/r + Z7) - uSl 

We plot the frequency dependence of the longitudinal conductivity in Fig [2l The collective 
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FIG. 2: The real imaginary and imaginary parts of axx, in units of cjq, for ^/uJc = 0.3 and u)cT = 3. 

cyclotron frequency ujc and the intrinsic, interaction induced damping frequency 7 were given 
in Eq. (ITTD . 
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For the thermoelectric response functions, we have 

{u + t/r) {aQ{fi/T)iij + t/r + 27) - tsp/ie + P)) 



X 



a^^(a;, k) 
a^yiuj, k) 



{u + i/r + i'^) — ujI 

ik (cj + z/r) (a; + i/r + i7 + ia;^/7) 

i — zvra — ctq 2 

^ {oj + i/T + — ojI 

[uj + i/t) {suJc/B + iaQ{j2/T){uj + i/T + 
{u + i/T + i^f' - ujI 

ik {uo + i/ r) (cj + i/r + 27 + iujl/ 7) 

1 — ivra — ctq 2 

^ {uj + i/t + i'j) — ujI 

axx{u,k) + 0{k^), 

s ujI + 7^ - 27 (g; + i/r) [1 - yip/jsT)] 

B (^LJ + i/r + 27)^ — o;^ 

ik {uo + ij r) (a; + i/r + 27 + iuJlj 7) 

X 1 — zTra — ctq 2 

(a; + i/r + 27) — u;^ 



(5.8) 

(5.9) 
(5.10) 



+ 0(fc' 



aa;j,(cj,fc) = dLxy{uJ,k) -^na 



kBaj 



{uj + i/r) 



+ 0{k^). 



(5.11) 



7" (cij + i/r + i7)^ — cUc^ 
Notice that the Onsager reciprocity a = a only holds for the response to homogeneous 
perturbations (i.e., /c = 0), while there are deviations at finite k. 
Finally, for the thermal conductivities, we have 



,P+e 



i^xxi^: k) 



{uj + i/r) + ctq^ {uj + i/r) (cj + i/r + i7) 



(u; + i/r + i7) -00'^ 
ik {lu + i/ t) {uj + i/r + i'y + iujl/ 7) 

X i — ZTTO; (Tq 2 

{uj + i/r + i'^) — ujI 

k 



i^xxi^: k) 

l^xyi}-^ 1 k) 



^ P + e uj + i/r 

-Zna-aQ —2 

J [00 + t/T + t'y) 

.(e + P) {u + i/T + iujlh) 



oot 



Oik'), 



T (cj + i/r + icj2/ -y)2 _ 



0(A;2 



(5.12) 
(5.13) 



B 



> 



^ - Ato-Q (7 



P+e 



- 2i 



P+e 



(tu + i/r + i7)) 



(cj + i/r + i7) - tUc^ 



ik (io + il r) (cj + i/r + i7 + itu^/ 7) 

1 — ivra — ctq 2 

{uo + i/r + i'y) — uo^ 



uo 



(u; + i/r + i7) -uoc 



^xyi^j k) 



(e + P) 



UOr. 



T {uj + i/r + iuol/'^Y - uo^ 



(5.14) 
(5.15) 



Note that the interaction-induced damping frequency for k, is uo'^/'j, and not 7. The former 
damping frequency also appears [ll| in the response functions for p^^ and be 
easily checked by inverting Eqs. (15.61) and (15.71) . 
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As in (15. 4p . in the presence of a magnetic field tlie tliermal conductivity in absence of 
currents is again independent of Coulomb interactions to lowest order in k/uj. 



C. Nernst effect 

An important thermo-electric response is the Nernst effect which measures the transverse 
electric field Ey that is established as a consequence of an applied longitudinal thermal 
gradient Va;T, in the absence of electrical currents. The ratio = Ey/ {—V xT) is called 
Nernst signal and is easily obtained from the coefficients determined above as cat = {a~^Q)xy 
It vanishes in the absence of a magnetic field and grows linearly for small fields B. The Nernst 
effect has become a popular measurement to characterize non-standard metals, such as in 



bismuth 26|] (where it was originally discovered 25|]), in other semimetals 27l |. in materials 



close to quantum critical points [28|], as well as in superconductors [29i]. All these systems 
share with undoped graphene the property of being far from a simple Fermi liquid. In the 
latter the so-called Sondheimer cancelation suppresses the Nernst signal, while it becomes 
very large in the systems mentioned above. By far the strongest Nernst signals (on the order 
of ImV/K for fields of IT) have been observed in bismuth whose band structure exhibits 
close similarities with graphene with which it shares the presence of nearly massless Dirac 
fermions. 

From the above formalism one easily obtains the full expression for the Nernst signal (at 
= 0) 

e kBTpico'J^ + l/rr + cor ^ ' 

where ks/e = 86.17/iV/K is its natural quantum unit. For small doping p, such that 
<^ min(7/co'c, 1), this result simplifies to 

In relatively clean samples close to quantum criticality this may exceed the quantum unit 
without violating the conditions for the applicability of hydrodynamics. 

In clean samples with uJcT ^ 1, and in the limit of large fields Uc <^ 7, one obtains the 
result 

kse + P 1 _ ks $g+p 
^ e ksTpu^T e {p/p,^Y{TT/h){B/B*y ^ ' ^ 

which decays inversely proportional to the field strength. 



VI. CONCLUSIONS 

Our main experimental predictions for the hydrodynamic cyclotron resonance in graphene 
are given by the frequency- dependent conductivities in Eqs. (15.61) and (15. 7p . with the fre- 
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quencies Uc and 7 specified in Eqs. (13. 5p and (13.81) . and the dimensionless function of the 
density $£+p estimated in Fig. [TJ This resonance occurs in a regime of weak magnetic fields 
where the Landau levels are not yet formed, and the dynamics is dominated by inelastic 
electron-electron collisions which occur at a rate ~ /i/fc^T. The electronic dynamics is 
"quantum critical", and the observation of such a resonance will offer a valuable oppor- 
tunity to explore quantum criticality. As has been argued elsewhere similar physics 
applies to a variety of systems in the vicinity of a quantum phase transition, including the 
superconductor-insulator transition in the cuprate superconductors. 
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